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A method to separate a Slater determinant wave function with a two-center neck struc- 
£\j | ture into spatially localized subsystems is proposed, and its potential applications are pre- 

sented. An orthonormal set of spatially localized single-particle wave functions is obtained by 
diagonalizing the coordinate operator for the major axis of a necked system. Using the local- 
1 ized single-particle wave functions, the wave function of each subsystem is defined. Therefore, 

defined subsystem wave functions are used to obtain density distributions, mass centers, and 
energies of subsystems. The present method is applied to separations of Margenau-Brink 
cluster wave functions of a + a, 16 + 16 0, and q+ 16 into their subsystems, and also to 
separations of antisymmetrized molecular dynamics wave functions of 10 Be into a + 6 He 
subsystems. The method is simple and applicable to the separation of general Slater deter- 
minant wave functions that have neck structures into subsystem wave functions. 

§1. Introduction 
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q ■ Nuclear systems often exhibit multi-center structures having spatially localized 

subsystems, as seen in phenomena such as cluster structures and fusion/fission pro- 
cesses. In cluster structures, each system consists of more than one spatially localized 
subsystems, called clusters, and intercluster motion is a degree of freedom used to 
describe the structure of a total system. In such multi-center systems, constituent 
nucleons of each cluster and the internuclear distance are useful parameters to specify 
cluster features. Identification of subsystems is also important for nuclear reactions, 
in which colliding nuclei are subsystems. In order to have an in-depth understand- 
ing of nuclear systems with subsystems, wave functions of each subsystem should be 
defined microscopically. 

For microscopic wave functions of quantum fermion systems, the Slater deter- 
minant description is often adopted as used in the Hartree-Fock (HF), the HF- 
Bogoliubov (HFB), and the time-dependent HF approaches, the Brink-Bloch model, 
and the antisymmetrized molecular dynamics (AMD) methodP'^ 1 '® Many applica- 
tions to studies of nuclear structures and reactions show that two-center structures 
have localized subsystems. However, a method to microscopically separate a Slater 
determinant into subsystem wave functions is not obvious, particularly when sub- 
systems overlap with each other. One of the popular ways to separate the density 
of a total system into those of subsystems is to employ the method of sharp-cut 
density separated by a boundary plane. However, the sharp-cut density is based on 
a classical picture and it is a nonphysical entity in quantum systems because density 
distributions of subsystems have singularities at the boundary plane. Recently, a 
method to localize quasi-particle states by unitary transformation applied to single- 
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particle orbits so as to maximize difference of occupations of left and right of the neck 
position in the HFB method (YG method).® In the YG method, localized single- 
particle wave functions depend on the neck position, which has ambiguities. In the 
cases of Margenau-Brink (MB) wave functions^ which are used in the Brink-Bloch 
model, and AMD wave functions, subsystem wave functions can be broadly defined 
as cluster wave functions or groups of single-particle Gaussian wave packets. How- 
ever, the subsystem wave functions in the overlapping region of two clusters contain 
nonphysical components that vanish because of the Pauli blocking effect, that is, 
antisymmetrization of nucleons in the total wave function. Saraceno et al. proposed 
canonical variables of the distance and relative momentum of the mass centers of 
the two clusters in an a— a system including the antisymmetrization effect between 
the a clusters (S method). However, wave functions of subsystems are not defined in 
the method,® and it is not applicable to systems that have one or more open-shell 
clusters. 

The aim of this paper is to propose a method applicable for any systems to 
define wave functions of spatially localized subsystems in a Slater determinant wave 
function by a linear transformation of single-particle orbits using an operator with 
respect to the internuclear relative coordinate. With this method, spatially localized 
single-particle wave functions are defined with no assumption and a total system 
is microscopically separated into subsystems. Applications to N = Z nuclei, 8 Be 
(a + a), 32 S ( 16 + 16 0), 20 Ne (a+ 16 0), and N / Z nucleus, 10 Be (a+ 6 He) are 
shown. We compare the internuclear distances defined by the present method with 
those defined by sharp-cut density, single-particle Gaussian centroids of AMD wave 
functions, and the S method. Energies of subsystems are also discussed. 

In £J21 we propose a method to define wave functions of subsystems. In £j3l 
the method is applied to MB wave functions for N = Z systems and AMD wave 
functions of 10 Be, and the density distributions of subsystems, internuclear distance, 
and energies of subsystems are shown. In §H characteristics of the method are 
discussed. Finally, conclusions are given in £}5j 

§2. Framework 

Suppose a Slater determinant wave function |$) having a two-center structure 

as 

\$)=A\<p 1) <p 2 , ,<pa), (2-1) 

where A is the antisymmetrization operator, A is the mass number of a total system, 
and cpx, a are single-particle wave functions. Here, the z axis is chosen to be the 
major axis of the total system. Note that the total wave function |$) is invariant 
except for the normalization under any linear transformations of single-particle wave 
functions {(f>{\ — > {<Pi} that give nonzero determinants. 

In order to separate the total system into two subsystems I and II, subsystem 
wave functions are defined as follows: First, single-particle wave functions are trans- 
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formed to diagonalize norm and coordinate operators z, 

=^2\<fj)cji, (2-2) 
j 

(ipi\ipj) = 5ij, (2-3) 

(<Pi\z\(pj) = ZiSij, (2-4) 

where z\ < Z2 < • • • < za- The eigenvalue Zi indicates the mean position of the i-th 
nucleon. If a single-particle wave function ipi is localized at a certain position, z; t 
can be regarded as the nucleon position in a semi-classical picture. Therefore, it is 
natural to classify <pi into two groups {<pi, <fi2, ■ ■ ■ , PAi} and {<pAi+i, i PA 1 +2, ■ • • , Pa} 
according to the distribution of z%. Then, wave functions and |$n) of the 
subsystems I and II, respectively, are described as 

=A\<pi,(p2,--- (2-5) 

|*n) = .A|¥>Ai+i, ¥>Ai+2, ■ • • , ^a)- (2-6) 

Here, and ^2 = A — A\ are regarded as the number of nucleons constituting 
subsystems I and II, respectively. 

Once single-particle wave functions are separated, the internuclear distance and 
expectation values of any operators (such as the density distribution and energy for 
each subsystem) can be calculated. The mass centers Ri and Rn of subsystems I 
and II are obtained as 



1 i=l,-,Ai 

2 i=A!+l,- ,A 

respectively. Then, the internuclear distance R is defined as the distance between 
the mass centers, 

J R=|R n -Ri|. (2-9) 
Density distributions p\ (r) and pn (r) of subsystems I and II are calculated by 

Pl (r)= ( 2 - 10 ) 

i=l,- ,Ai 

pn(r)= ( 2 ' n ) 

i=A 1 +l,- ,A 

respectively. Here, pi(r) is the density distribution of the i-th nucleon, 

Pi (T) = (<p i \6(r-r)\<p i ). (2-12) 
Hamiltonian H n for the subsystem n is defined as, 

H n = f n + V n -f Gn , (2-13) 
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where T n , V n , and Tq u are the kinetic, potential, and mass-center motion terms, 
respectively, for the subsystem n. The kinetic term T^ = T n — Tcj, n with the center- 
of-mass correction is described as a two-body operator (see Appendix). Energy E n 
of the subsystems n is obtained as 

E n = {<S> n \H n \<f> n ). (2-14) 

In the present study, the Volkov No. 2 force is used as the effective interaction, and 
the Majorana parameter is set to 0.59 to adjust the ground-state energy of 16 0. 

§3. Results 

3.1. N = Z nuclei 

In this section, the current separation method is applied to the MB wave func- 
tions. In the MB wave function for a system consisting of two clusters, clusters I 
and II are expressed by the shell-model configurations with centers at (0,0, — -jr-d) 
and (0, 0, +-^-d), respectively. The parameter d specifies the degree of cluster devel- 
opment. MB wave functions for two-cluster systems, a + a, a + 16 O, and 16 O + 16 O, 
are considered, and subsystem density distributions, internuclear distances, and the 
energy of each cluster are calculated by the present method. 

3.1.1. 8 Be(a + a) 

In a MB wave function |3> Q+Q ) with an a + a cluster structure, the a clusters I 
and II are expressed by the (0s) shell model wave function shifted to the positions 
— d/2 and d/2, respectively, as 



\<S>a+a) 


= A\<fil, <P2,-" 1 &8), 


(3-1) 


\<Pl,s) 


= |01,2 ®pt)i 


(3-2) 


\<P2,6) 


= |01,2 ®Pl), 


(3-3) 


l^3j) 


= |01,2 (8> Tit). 


(3-4) 


l^4,8> 


= |0i,2 ® n |), 


(3-5) 


(r|0l,2> 




(3-6) 



where the d = (0, 0, d), p and n are a proton and neutron, respectively, the f and \ 
indicate up and down spins, respectively, and the v indicates the width parameter 
of Gaussian wave packets. In the large d limit, the antisymmetrization effect is 
small enough and the mean position of a clusters I and II are z Fd/2, respectively. 
Therefore, in a semi-classical picture, the parameter d is regarded as the internuclear 
distance. However, when the spatial parts 0i and 02 have a non-negligible overlap 
with each other, d does not represent the internuclear distance because the single- 
particle wave functions in the overlapping region are modified by the Pauli blocking 
effect. 

The transformation given by Eq. (|2-ip - (|2-4p obtains analytical forms of spatially 
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localized single-particle wave functions in the a + a system as 



1^1,5) 


= |01,2 <8> 


Pt), 


|V2,6> 


= 1^1,2 <8> 




1^3,7) 


= |01,2 ® 


nt) 


l^4,8> 


= |01,2 ® 


n |) 


<r]0i, 2 ) 


= c+e"^ 






1 /2zv 


3 / 
\ 4 / 



(3-7) 
(3-8) 
(3-9) 
(3-10) 




Eigen values {zj} of I are 



7T 



*1,2,3,4 = ■ (3-13) 

2V1 - e-^ d 

Z 5 ,6,7,8 = + , ^ = , (344) 

2V1 - e-^ d 



which are scaled by l/yl-e"^ when compared to the positions of centroids of 
wave packets =Fd/2. The internuclear distance R between a clusters is, therefore, 

R= , d = ■ (3-15) 

The internuclear distance R tends to z/ -1 / 2 in the d — > limit. 

Figure [T] shows the spatial parts of the spatially localized single-particle wave 
functions 0i 5 2 as well as Gaussian wave packets 01,2 at d = 4, 2, and 1 z^ -1 / 2 , and 
the d — > limit. In the case of d = 4 i/ -1 / 2 , 0j almost coincides with 0j, whereas 0j 
has an obvious node in the d < 2z^ -1 / 2 region where 0i and 02 overlap. The wave 
functions 0i and 02 are spatially localized even for d — > 0. 

Figure [2] shows linear density distributions p xy i,u(z) of a clusters I and II defined 
by the integrated subsystem densities with respect to x and y, 

Pxy\,u{z) = JJdxdy /9i,n(r). (3-16) 

Density distributions of four single-particle wave functions in each cluster are the 
same because of spin and isospin saturation. Unlike the sharp-cut separation, sub- 
system densities are always smooth. Density distributions of a clusters I and II are 
spatially localized, even in the case of d — > 0. 

Figure [3] shows internuclear distances R by the present method and R p obtained 
by the sharp-cut separation as functions of distance parameter d between the cen- 
troids of wave packets. R p is calculated analytically as 



Rn 



v 7 2ttv sinh (\vd 2 ) 



(3-17) 
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Fig. 1. Amplitudes of a 2-part of single-particle wave functions <j>i (solid) and <j>2 (dashed), where 
a z-part of <f>i is defined as <j>i/ , in a + a MB wave functions are plotted 

with the parameters d = 4 (upper left), 2 (upper right), and 1 u" 1 ^ 2 (lower left), while d->0 
(lower right) are plotted as functions of z. The Gaussian wave packets for the original single 
particle wave functions </>i,2z (dotted and dot-dashed, respectively) before antisymmetrization 



are also shown. Units of the amplitudes and z are in v 1 ^ and 



-1/2 



respectively. 



In d > 2 u- l l 2 , R and R p are similar and they are close to d. In d < 2 v x / 2 , where 
a clusters overlap with each other, R is smaller than R p . 

Figure [5] shows energies of the a clusters in the 8 Be as function of internuclear 
distance R. Width parameter v is set to 0.27 fm~ 2 . In R > 5 fm region, the a 
clusters are in the ground state, but in R < 5 fm region, the a clusters are distorted 
by the antisymmetrization effect between clusters, and their excitation energies go 
up smoothly. 

3.1.2. 20 Ne (a+ 16 0) 

Figure shows linear density distributions for the a and 16 O clusters, as well 
as spatially localized single-particle orbits of a and 16 clusters in a+ 16 MB wave 
functions with d = 8,4,2, and 1 fm. Width parameter v is set to 0.16 fm -2 . At 
the small d region, a and 16 clusters have nodes in the overlap region because 
of antisymmetrization effects between clusters, but they are still spatially localized. 
Subsystem densities have no singularity. 

Figure [6] shows the internuclear distance R and R p between a and 16 clusters 
defined by the present method and sharp-cut separation, respectively. R is smaller 
than R p in the d < 6 fm region, and they are similar in the d > 6 f m region. 



Separation into subsystems 



7 






4.0 r 




3.5 - 




3.0 - 


>i 2.5 - 




2.0 - 


<o 


1.5 - 




1.0 - 




0.5 - 



0.0 







single I 


-\d = 2 v- 1 ' 2 




single TT 






_ /'i 






/ \ .'ii 






/ \ total 




\ / 





1 
.y-V2l 



1.0 
3.5 
3.0 - 



a 2 -° 

S 

XJ 1.5 

| 1.0 

a 

W 0.5 



(1.0 



-2 -1 

Z If 



1 

1/21 





A single I 


- 1 d — > | / \ 


£ \ single II - 




1 \ Pi 




1 Pll 




1 total 


/ ! 









Fig. 2. Linear density distributions of a clusters I and II, spatially localized single-particle wave 
functions, and the total system in a + a MB wave functions with the parameter d = 4 (upper 
left), 2 (upper right), 1 u" 1 ^ 2 (lower left), and d — > (lower right) are plotted as functions of z. 
Thin dot-dashed and dotted lines represent single-particle wave functions for a clusters I and 
II, respectively. Thick dot-dashed and dotted lines represent a clusters I and II, respectively. 
Solid lines represent total wave functions. Units of linear densities and z are in v 1 ^ 2 and v~ x ^ 2 , 
respectively. 



Figure [7] shows energies of the a and 16 clusters in the 20 Ne as functions of 
internuclear distance R. In R > 6 fm region, the a and 16 O clusters are in the 
ground state, but in R < 6 fm region, the a and 16 clusters are distorted by the 
antisymmetrization effect between clusters, and their excitation energies go up. 

3.1.3. 32 S( 16 0+ 16 0) 

Figure [8] shows linear density distributions of 16 O clusters and spatially localized 
single-particle wave functions in 16 + 16 MB wave functions with d = 10, 6, 5, 
and 2 fm. Width parameters v are set to 0.16 fm -2 . Each total density distribution 
shows a neck between two clusters, and each single-particle wave function is well 
localized. Therefore, nucleons are separated into two 16 clusters. The densities 
of the total system and the two 16 clusters are also shown in the figure. In the 
touching region such as d = 5 and 6 fm, the density profile for each 16 O cluster 
defined in the present method also shows reasonable shapes. Subsystem densities 
have no singularity, and are smoothly damped in the overlap region. 

The internuclear distance R determined by the present method is plotted as a 
function of d in Fig. The distance R p defined by the sharp-cut density is also 
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Fig. 3. Internuclear distance defined by 
the present method (solid) and sharp-cut 
separation (dashed) as functions of dis- 
tance between centroids of wave packets 
in units of v^ 1 ^ 2 . 




1 2 3 4 5 6 
internuclear distance fm 

Fig. 4. Energies of the a cluster in 8 Be are 
plotted as a function of internuclear dis- 
tance in the unit of MeV/A (solid), where 
A is the mass number of each cluster. 
The dashed line shows the ground-state 
energy of the a cluster. 




Fig. 5. Linear density distributions of a and 16 clusters and those for the total system in the 
a+ ie O MB wave functions with the parameter d = 8, 4, 2, and 1 fm. Thin dot-dashed and 
dotted lines represent single-particle wave functions for a and 16 clusters, respectively. Thick 
dot-dashed and dotted lines represent a and 16 O clusters, respectively. Solid lines represent 
total wave functions. Units of linear densities and z are in fm -1 and fm, respectively. 
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Fig. 6. Internuclear distance between a 
and 16 O clusters defined by the present 
method (solid) and sharp-cut separation 
(dashed) as functions of distance between 
centroids of wave packets. Units are in 
fm. 
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Fig 



7. Energies of the a (solid) and 16 
(dot-dashed) clusters in 20 Ne are plotted 
as functions of internuclear distance in 
the unit of MeV / A, where A is the mass 
number of the cluster. The dashed and 
dotted lines shows the ground-state en- 
ergy of the a and 16 clusters, respec- 
tively. 



shown for comparison. R equals to d in the large d region, and as d decreases, it 
deviates from d toward the R > d region. In the d = limit also, the internuclear 
distance R has a finite value because of the Pauli blocking effect. The behavior of R 
is qualitatively similar to that of R p . However, quantitatively, R is slightly smaller 
than R p because the left and right nuclei have an overlap density in the present 
definition instead of the sharp-cut density. 

Figure [TOl shows energies of the 16 O clusters in the 32 S as a function of internu- 
clear distance R. In R > 7 fm region, the 16 clusters are in the ground state, but 
in R < 7 fm region, the 16 O clusters are distorted by the antisymmetrization effect 
between clusters, and their excitation energies go up. 

3.2. N + Z: 10 Be (a+ e He) 

In this section, the present method is applied to AMD wave functions for N ^ 
Z nucleus 10 Be that is calculated by energy variations after parity and angular 
momentum projection. 61 In the AMD model, an intrinsic wave function is expressed 
by a Slater determinant of single-particle wave functions each of which is written by a 
Gaussian wave packet. Gaussian centers of all nucleons are treated as independently 
varying parameters, and therefore, no clusters are assumed. One of the advantages of 
the AMD method is that the formation and/or breaking of clusters can be described 
in the framework. Since clusters are not assumed in the AMD model, a grouping of 
nucleons into two clusters is not necessarily obvious, unlike the MB cluster model 
wave functions. 

Analyzing the mean position Zi of nucleons for intrinsic wave functions deter- 
mined by the AMD method, single-particle wave functions tpi can be separated into 
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Fig. 8. Linear density distributions of 16 clusters I and II, and the total system in 16 0+ 16 
MB wave functions with the parameter d = 10, 6, 5, and 2 fm are plotted as functions of z. 
Thick dot-dashed and dotted lines represent single-particle wave functions for 16 O clusters I and 
II, respectively. Thin dot-dashed and dotted lines represent 16 clusters I and II, respectively. 
Solid lines represent the total wave functions. Units of linear densities and z are in fm -1 and 
fm. 
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Fig. 9. Internuclear distance between 16 O 
clusters defined by the present method 
(solid) and sharp-cut separation (dashed) 
as functions of the distance between cen- 
troids of wave packets. Units are in fm. 
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8 9 10 Fig. 10. Energies of the 16 cluster in 32 S 
are plotted as a function of internuclear 
distance in the unit of McV/A (solid), 
where A is the mass number of the clus- 
ter. The dashed line shows the ground- 
state energy of the 16 cluster. 
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Fig. 11. Linear density distributions of intrinsic wave functions of dominant components of the 
J n = 0^ (left) and 0^ (right) states in 10 Be. Thin dot-dashed and dotted lines represent single- 
particle wave functions for a and 6 He clusters, respectively. Thick dot-dashed and dotted lines 
represent a and 6 He clusters, respectively. Solid lines represent total wave functions. Units of 
linear densities and z are in fm -1 and fm. 



two subsystems I and II. Four nucleons have Zi in the left of the neck position and 
form an a cluster (subsystem I). Six nucleons having Z{ in the right of the neck can 
be classified in subsystem II corresponding to a 6 He cluster. Figure [TT1 shows linear 
density distributions of intrinsic states of dominant components of the J n = Of and 
states in 10 Be, which are separated into a and 6 He clusters. In both states, total 
linear density distributions have necks at z ~ fm, and single-particle orbits localize 
left and right parts. Therefore, the entire 10 Be system can be separated into two 
clusters using the present separation method. 

§4. Discussions 

The previous section shows that the present separation method involving the 
diagonalization of a spatial operator of the major axis works well for the separation 
of a Slater determinant wave function into spatially localized wave functions of sub- 
systems. The results of applications to MB wave functions with a + a, a + 16 0, 
and 16 O + 16 O structures are shown in Figs. [21 [5J and El respectively. With this 
method, we obtain subsystem wave functions with no singularity. Because of non- 
singular wave functions, energies of subsystems do not diverge as shown in Figs. HI 
[71 and [TUl This is one of the advantages over the sharp-cut method. Furthermore, 
the present method provides single-particle wave functions of subsystems, which are 
useful in the microscopic analysis of structures of subsystems, and is applicable to 
studies of structural changes of subsystems in cluster structures and nuclear reac- 
tions, depending on the internuclear distance. In the applications involving AMD 
wave functions of J 77 = Of and states in 10 Be, although the existence of clusters 
is not assumed, the proposed method is proven to be useful in the separation of 10 Be 
wave functions into subsystems a and 6 He, when the neck positions and the single- 
particle orbits (Fig. [11]) are compared. The proposed method is applicable to the 
separation of Slater determinant wave functions with neck structures into spatially 
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Fig. 12. Internuclear distance for the a + a MB wave functions specified by the parameter d. The 
distance R defined in the present method, Rs defined in the S method, and R p defined by the 
sharp-cut density are plotted. Units are in v -1 ' 2 . 



localized subsystems. 

As mentioned above, the present method gives wave functions of subsystems 
with which expectation values of operators can be calculated. For example, the 
calculation of mass centers of subsystems determines the internuclear distances of a 
+ a, a + 16 0, and 16 + 16 systems, as shown in Figs. El El andO respectively. The 
present method is useful to study cluster structures and nuclear reactions because 
the internuclear distance is an important degree of freedom in these phenomena. 

As an internuclear distance including the antisymmetrization effects between 
clusters, Saraceno et al. proposed the internuclear distance as one of canonical 
variables for an a + a system.'' ' 



Rs = ^ Y ±^d=VTT^R. (4-1) 

We compare Rs with R and discuss properties of these internuclear distance that 
include the antisymmetrization effects between clusters. Although definitions of 
internuclear distances R and Rs are different, the form of R in the present definition 
resembles that of Rs in the denominator, which is related to the antisymmetrization 
effect, but is smaller than Rg because of the factor yl + e~ ud2 (Fig. [T2"|) . and in the 
d = limit, R = Rs/V2- In the present method, single-particle orbits are defined to 
be completely orthogonal. Under this condition, the internuclear distance is smaller 
than Rs- With the increase in d, both the distances R and Rs approach d. 

Advantages of the present method are as follows. First, subsystem wave func- 
tions can be obtained. They are expressed by an orthogonal set of single-particle 
orbits so that the Fermion feature of constituent nucleons is respected. The internu- 
clear distance given by our method is the largest value under this condition. With 
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the wave functions, subsystem properties such as the subsystem energy can be dis- 
cussed. This is an advantage over the S method and the sharp-cut density method. 
Second, our method is applicable to any systems even in cases of proton- or neutron- 
excess systems or systems having one or more open-shell clusters or both. Third, 
our proposed method obtains spatially localized single-particle wave functions with 
no ambiguity in a simple linear transformation, which can be an advantage over the 
YG method. In the YG method, the neck position has to be determined to obtain 
localized single-particle wave functions. The determination of the neck position may 
have less ambiguity in a case of a well separated system at a scission point, but it 
may have ambiguities in particular when two nuclei overlap to each other. Moreover, 
in a case of much proton- or neutron-rich systems, the neck position is unclear and 
neck positions of protons and neutrons can be different. 

§5. Conclusions 

In this paper, we proposed a method to separate a Slater determinant wave 
function with a neck structure into spatially localized subsystems. Our method is 
applied to the MB wave functions of a + a, a + 16 0, and 16 + 16 systems, and 
the AMD wave functions of dominant components of the J 71 " = Of and 0^ states in 
10 Be. 

The proposed method obtains wave functions of subsystems with no singularity, 
which is an advantage over the sharp-cut method. Using the obtained subsystem 
wave functions, we calculated the expectation values of several operators for subsys- 
tems. For example, the internuclear distance is well defined by the calculated mass 
centers of the subsystems, and energies of subsystems do not diverge. 

The proposed method is simple and applicable to a wide variety of approaches 
that are based on the Slater determinant wave function, e.g., the HF method. It is 
useful for the analysis of systems that have spatially localized subsystems in phe- 
nomena such as cluster structures and nuclear reactions. 
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Appendix 

The kinetic term T' = T — Tq with the center-of-mass correction is described as 
a two-body operator as follows, 

f = f - f G 

= y- p[_ _ (Ei Pi) 2 
^ 2m 2Am 

i 
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2m EP* 2 2Am^ t>i ' i>j 

% ij 



;2 , i2 



1 ^ Pt + Pj , A 

2^ 9 - 1* • Pi 



2Am ^ \ 2 



(Pi - Pi 



= —Y 

v 

= V%#. (A-l) 



By a similar procedure, the mean squared radius r 2 with the center-of-mass 
correction is described as 



J / % <3 
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